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ABSTRACT 

The decomposition of a signal on the sphere with the steerable wavelet constructed 
from the second Gaussian derivative gives access to the orientation, signed-intensity, 
and elongation of the signal's local features. In the present work, the non-Gaussianity 
of the WMAP temperature data of the cosmic microwave background (CMB) is an- 
alyzed in terms of the first four moments of the statistically isotropic random fields 
associated with these local morphological measures, at wavelet scales corresponding 
to angular sizes between 27.5' and 30° on the celestial sphere. While no detection is 
made neither in the orientation analysis nor in the elongation analysis, a strong de- 
tection is made in the excess kurtosis of the signed-intensity of the WMAP data. The 
non-Gaussianity is observed with a significance level below 0.5% at a wavelet scale 
corresponding to an angular size around 10°, and confirmed at neighbour scales. This 
supports a previous detection of an excess of kurtosis in t he wavelet coefficie nt of the 
WMAP data with the axisymmetric Mexican hat wavelet (jVielva et al.ll200^ . Instru- 
mental noise and foreground emissions are not likely to be at the origin of the excess 
of kurtosis. Large-scale modulations of the CMB related to some unknown system- 
atics are rejected as possible origins of the detection. The observed non-Gaussianity 
may therefore probably be imputed to the CMB itself, thereby questioning the basic 
inflationary scenario upon which the present concordance cosmological model relies. 
Taking the CMB temperature angular power spectrum of the concordance cosmologi- 
cal model at face value, further analysis also suggests that this non-Gaussianity is not 
confined to the directions on the celestial sphere with an anomalous signed-intensity. 

Key words: methods: data analysis, techniques: image processing, cosmology: ob- 
servations, cosmic microwave background 



1 INTRODUCTION 

Among other cosmological observations, the recent data 
of the cosmic microwave background, in particular those 
released by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite experiment, have played a central role 
in defining a concordance model highlighting a flat ACDM 
Universe with a primordial phase of inflation ([Bennett et al 



2003 : ISpergel et al ] |2003l : lHinsh"aw et a ] |2007l : ISpergel et al 
20071). In that framework, the flat Universe is flUed in with 



cold dark matter (CDM) and dark energy in the form of 
a cosmological constant (A), in addition to the standard 
baryonic and electromagnetic components. The cosmological 
principle is assumed, which postulates global homogeneity 
and isotropy. In the basic inflationary scenario considered, 
the large scale structure and the cosmic microwave back- 



ground fluctuations are also assumed to arise from Gaus- 
sian quantum energy density perturbations around the ho- 
mogeneous and isotropic background during the inflationary 
phase. The cosmological parameters of the concordance cos- 
mological model are determined with a precision of the order 
of several percent, but the hypotheses on which the model 
relies need to be thoroughly challenged. 

In that context, the CMB constitutes a realization of a 
random field on the sphere. The cosmological principle and 
the basic infiationary scenario respectively imply the sta- 
tistical isotropy and the Gaussianity of this random field. 
These basic hypotheses of the concordance cosmological 
model may consequently be questioned through the CMB 
analysis. 

The statistical isotropy of the CMB tempera- 
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ture field has already been largely questioned in 
the analysis of the WMAP data. Firstly, a North- 



tected (Eriksen et al. 2004allbl. 12005; Hansen et al.ll2004a bl: 


Donoehue & Donoghue 


2005: 


Land & Maeuciio l2005al: 


Bernui et al.l 2009. 1200 


71: iSoergel et al.. i2007i: .Eriksen et al. 


20071: iMontesen'n et all |2007|). 


Secondly, an anomalous 



aliRnme nt of the lowest multipoles of the data was ob- 



CoDi et al. l2004l: Katz & Weeks 


2004: Bielewicz et all 
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: iFreeman et al. 20061: 


Abramo et al. 20061). Finally, wavelet analyses have also re- 



ported statistical isotropy anomalies related to the signed- 
intensity of local CMB features, as well as anomalies re- 
lated to the alignment of local CMB features toward spe- 
cific direc tions on the celestial sphere (| Wiaux et al.ll2006al : 
IVielva et al. 2006, 2007). 

The Gaussianity of the CMB temperature field has 
also been largely questioned in the analysis of the WMAP 
data. Firstly, departures fr om Gaussianity were detected 
using statisti c s of e xtrema (|Lars on fc Wandclt' 2004, 2005'; 
iToieiro et al.l |2006| ) , bispectra ^tLM!A , ^iJ^^^}^^}\9. jOOSg); 



phase correlations (I Chiang et al 



l2006l : IColes et al 



2004 



Naselskv et al 



Chiang fc Naselskvl 
20051 ). Minkowski 



function als l|Parkll2004l: lEriksen et al.ll2004bl ). and local cur 



vature (|Hansen et al.l 12004a''). Secondly, wavelet analyses 
have also reported non-Gaussian deviations. An excess of 
kurtosis in the wavelet coefficient of the WMAP tempera- 
ture data with the axisymmetric Mexican hat wavelet on 
the sphere was found at wavelet scales corresponding to an- 
gular sizes on the celestial sphere aro und 10°, and localized 
in the southern galactic hemisphere (jVielva et al.|[2004i ). A 
cold spot (i.e. with negative wavelet coefficients) was iden- 
tified at (0,^) = (147°, 209°), with 9 € [0,7r] and ip e [0, 27r) 
respectively standing for the co-latitude and longitude in 
galactic spherical coordinates, and considered to be a good 
candidate to explain the observed deviation. The confirma- 
tion that the cold spot is anomalous was provided, still with 
the axisymmetric Mexican hat wavelet, in terms of its area 
dCruz ct al. 2005). The detection was further confirmed with 



various wavelets and various statistics (IMukheriee fc Wang 


2004 


: Cavon et al. 20051: 


McEwen et al. 20051: Cruz et al. 


2006 


.12007: McEwen et al. 


2006f). Notice that the cold soot 



identified also certainly represents a departure from statisti- 
cal isotropy, in terms of a North-South asymmetry in galac- 
tic coordinates. 

By essence, wavelet analyses present the particular ad- 
vantage of probing, not only the scale but also the localiza- 
tion of the features consti tuting the CMB on the celestial 
sphere (|Wiaux et al.ll2005l ). Steerable wavelets also provide 
morphological measures of the local f eatures, such as orien - 
tation, signed-intensity, or el ongation (iMcEwen et al.ll2007t ). 
at a low computational cost (W iaux et al .l2006bl ). They were 
used to probe the statistical isotropy of the WMAP CMB 
temperature data in the previously quoted signed-intensity 
and alignment analyses. They were also used to probe the 
Integrated Sachs-Wolfe effect through the correlation of 
WMA P data and large scale structure data (McEwen ct aU 
I2OO7I ). In the present work, a further insight into the CMB 
temperature non-Gaussianity is provided through a steer- 
able wavelet analysis of the WMAP data. 

In Section (2] we present the methodology adopted. In 



Section[3l we present the results of the WMAP data analysis. 
In Section [4] we study systematic effects as a possible origin 
of the detections. In Section [S] we discuss the origin of our 
detection and its detailed interpretation. We finally conclude 
in Section [6l 



2 METHODOLOGY 

In this section, we firstly recall the formalism for the anal- 
ysis of signals on the sphere with steerable wavelets, as well 
as the local morphological measures of orientation, signed- 
intensity, and elongation, defined from the steerable wavelet 
constructed from the second Gaussian derivative. Secondly, 
we explicitly describe the statistics for the non-Gaussianity 
analysis on the random fields associated with the local mor- 
phological measures of the CMB temperature field. These 
statistics are simply the first four moments of the random 
fields considered. 



2.1 Steerable wavelets and morphological 
measures 

We consider the three-dimensional Cartesian coordinate sys- 
tem (o, ox, oy, oz) centered on the unit sphere, and where the 
direction oz identifies the North pole. Any point ua on the 
sphere is identified by its corresponding spherical coordi- 
nates {9, tp), where 9 € [0, tt] stands for the co-latitude, and 
ip G [0, 27r) for the longitude. 

Firstly, we briefly summ arize the formalism of steerable 
wavelets on the sphere S'^ (I Wiaux et al.|[2005l ). Any fllter 
invariant under rotation around itself is said to be axisym- 
metric. By definition, any non-axisymmetric, or directional, 
filter is steerable if a rotation by x € [0, 2ii) around itself 
may be expressed in terms of a finite linear combination of 
M non-rotated basis filters ^m'- 



H = (X) *r, 



(1) 



where the weights km{x), with 1 ^ m ^ M, are called 
interpolation functions. The analysis of a signal F with a 
given wavelet ^ simply defines a set of wavelet coefficients 
{ijJo,x,a), which result from the directional correlation 
between F and the wavelet dilated at any scale a, v&a. In 
other words these wavelet coefficients are defined by the 
scalar product between the signal and the wavelet dilated 
at scale a, rotated around itself by x, and translated at any 
point ujQ on the sphere, also denoted ^uja,x,a'- 

Wi{cuo,X,a) = {^^„,x.a\F)= f dQKo,x,aHFH. (2) 

The * denotes complex conjugation. The wavelet coefficients 
of a signal therefore characterize the signal at each scale a, 
orientation x, and position uo- 

In the present work, we consider the second Gaussian 
derivative wavelet (2GD), ■^^a'^sau) ^ which is obtained by a 
stereographic projection of the second derivative in direction 
i of a Gaussian in the tangent plane at the North pole. The 
filter obtained is a steerable wavelet on the sphere which 
may be rotated in terms of three basis filters (M = 3): the 
second derivative in the direction x itself, the second deriva- 
tive in the direction y, and the cross-derivative. Notice that 
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the value of the scale a identifies with the dispersion of the 
Gaussian in units of 2tan(S/2). The angular size of the 2GD 
is defined as twice the half- width of the wavelet, where the 
half-width is defined by 6hw = 2 arctan(a/2), which is closely 
approximated by a at small scales. 

Secondly, we recall that the 2GD gives access to three lo- 
cal morphologic al measures of orient ation, signed-intensity, 
and elongation (jMcEwen et al.| [2007l). By linearity, the re- 
lation of steerability ([T]) is automatically transferred on the 
wavelet coefficients of F. Consequently, at each scale a and 
at each position ljo, the orientation xo(i^o, a) that maximizes 
the absolute value of the wavelet coefficient, can easily be 
computed, with an infinite theoretical precision. It corre- 
sponds to the local orientation at which the wavelet best 
matches the local feature of the signal. As the 2GD is in- 
variant under rotation around itself by tt, orientations may 
arbitrarily be constrained in a range of length tt, and as the 
2GD oscillates in the tangent direction x, it actually detects 
features aligned along the tangent direction y. The local ori- 
entation of the feature itself, {ijjQ,a), is therefore defined 
in terms of xo = Xo('^o,a) as: 



2 



{too, a) 



TT 

Xo+2 



(3) 



The wavelet coefficient itself at scale a, position ujo, and in 
direction xo , defines to so-called signed-intensity of the local 
feature: 

l'' {uo,a) = W^g2{uo,xo,a) . (4) 
The elongation of local features is explicitly defined by 



0^ E^{uo,a) = l 



(a;o,Xo + f,a) 



<1. (5) 



Numerical tests performed on elliptical Gaussian-profile 
features show that this elongation measure increases 
monotonously in the range [0, 1] with the intrinsic eccen- 
tricity e £ [0, 1] of the features. While it is possible to define 
alternative elongation measures, these numerical tests also 
indicate that the chosen definition is not an arbitrary mea- 
sure of the non-axisymmetry of local features, but represents 
a rough estimate of the eccentricity of a Gaussian-profile lo- 
cal feature. 

In summary, the analysis of signals with steerable 
wavelets is interesting in several respects. Firstly, the wavelet 
decomposition enables one to identify the scales a of the 
physical processes which define the local feature of the sig- 
nal at each point ujo. Secondly, the steerability theoretically 
gives access to local morphological measures. For the 2GD, 
the orientation, signed-intensity and elongation of local fea- 
tures are defined. Finally, from the computational point of 
view, the calculation of a directional correlation at each anal- 
ysis scale is an extremely demanding task. The relation of 
steerability is essential to reduce the complexity of calcula- 
tion of the wavelet coefficients when local orientations are 
considered l|Wiaux et al.ll2006bl ). 



2.2 Statistics for non-Gaussianity 

In the context of the concordance cosmological model, the 
CMB temperature represents a realization of a statistically 
isotropic and Gaussian random field on the sphere. The 



WMAP data are also contaminated by noise and foreground 
emissions. The statistical analysis is performed by compar- 
ison of the data with simulations. The noise present in the 
data is simulated and the regions of the sky in which the 
data are too much contaminated by foreground emissions 
are masked, and excluded from the analysis. Typically, a 
non-Gaussianity analysis is performed through the evalua- 
tion of global estimators computed as simple averages on the 
whole part of the celestial sphere where the data are consid- 
ered to be valid, explicitly assuming the statistical isotropy 
in the corresponding part of the sky. Any anomaly between 
the data and the simulations is consequently interpreted as 
a departure of the data from Gaussianity. 

We consider the statistically isotropic real-valued ran- 
dom fields on the sphere associated with the local morpho- 
logical measures of orientation, signed-intensity, and elon- 
gation of the CMB temperature field T at each wavelet 
scale a: X(cjo,a), with X = {D'^ j"^ , E"^}. The statis- 
tics estimated for the subsequent non-Gaussianity analysis 
are simply moments of the first four orders. The first two 
are the mean fi'^ (a) — (X(cjo,a)), and variance ^"^(a) = 
{[X{ujo,a) — jj,'^ (a)]'^)^^^ . The third-order moment is the 
skewness (a) = {[X{Luo,a) - fi^ (a)]^) /[a^ (a)]^ . The 
skewness measures the asymmetry of the probability den- 
sity function, and hence a deviation relative to a Gaussian 
distribution. Positive and negative skewnesses are respec- 
tively associated with larger right and left distribution tails. 
The fourth-order moment considered is the excess kurtosis 
A'-^(a) = {[X{wo,a) - fi^ {a)]^) /[a^ (a)]* - 3. The kurtosis 
measures the peakedness of the probability density function 
relative to a Gaussian distribution. Positive and negative ex- 
cess kurtoses are respectively associated with distributions 
more and less peaked than a Gaussian distribution. These 
four moments are independent of the point uo because of the 
statistical isotropy. The corresponding estimators computed 
by averages over the sphere are 



(a) 



K"" (a) 



Z — 1 

[X(c.«,a)-M^ (a) 



1/2 



- -y 



(a) 

X (a) 



■5^ (a) 



(6) 



At each wavelet scale a, Na stands for the total number 
of valid pixels outside a given exclusion mask Ma which, 
by definition, identifies the pixels to be excluded from the 
analysis (see Subsection 13. ip . The values cJq'^ identify the 
center of these valid pixels. 

Let us emphasize here that, even under the assump- 
tion of Gaussianity of the CMB temperature field, postu- 
lated for the simulations, none of the random fields associ- 
ated with these local morphological measures at each scale 
is intrinsically Gaussian. This is simply due to the non- 
linearity of the definitions ([3| for the orientation, Q for the 
signed-intensity, and ((5)l for the elongation. In particular. 
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by statistical isotropy, the measure of orientation {loq, a) 
should be uniformly distributed at each point ljq and at each 
wavelet scale a. 

For each local morphological measure, statistics, and 
wavelet scale, the value obtained for the data can be com- 
pared to the corresponding values for the simulations. The 
percentiles corresponding to specific cumulative probabili- 
ties p in the simulations considered are calculated for a first 
comparison with the value of the data. The percentile as- 
sociated with p = 50% defines the median value. Cumula- 
tive probabilities p = {15.865%, 84.135%} are considered, 
which formally correspond to the percentiles at one stan- 
dard deviation (la) from the mean in a Gaussian distribu- 
tion. They define a first, innermost, region for the distri- 
bution of percentiles around the median value. The exact 
values considered reflect the maximum precision allowed by 
our sample of ten thousand simulations. Cumulative prob- 
abilities p — {2.275%, 97.725%} are also considered, which 
formally correspond to the percentiles at two standard devi- 
ations (2cr) from the mean in a Gaussian distribution. They 
define a second, middle, region for the distribution of per- 
centiles around the median value. Again, the exact values 
considered reflect the maximum precision allowed by our 
sample of ten thousand simulations. Cumulative probabil- 
ities p — {0.5%, 99.5%} are finally considered, defining a 
third, outermost, region for the distribution of percentiles 
around the median value. If the value of the data for a statis- 
tics at a given scale is higher (lower) than the median value 
obtained from the simulations, the significance level of a de- 
tection is simply defined as the fraction of simulations with 
a higher (lower) value than the data. The lower the sig- 
nificance level, the stronger the detection. Typically, in the 
following, a significance level below 0.5%, corresponding to 
values outside the outermost region for the distribution of 
percentiles around the median value, will be associated with 
a strong detection. 



3 WMAP ANALYSIS 

In this section, we firstly describe the pre-processing proce- 
dure applied to the three-year WMAP CMB temperature 
data, as well as to the corresponding simulations produced 
from the concordance cosmological model. Secondly, we ex- 
pose the results of the application of the non-Gaussianity 
analysis defined in the previous section on the three-year 
WMAP co-added CMB data, notably highlighting a strong 
detection in the excess kurtosis of the signed-intensity. 



3.1 Data and simulations 

Firstly, the following pre-processing procedure is applied to 
the three-year WMAP CMB temperature data before the 
non-Gaussianity analysis. The original maps of the eight 
WMAP radiometers at the Q, V, and W frequencies (Ql and 
Q2 at 41 GHz, VI and V2 at 61 GHz, and Wl, W2, W3, and 
W4 at 94 GHz) are corrected for foreground emissions con- 
tamin ation by a template fitting technique (jSpergel et al.l 
[2003). The resulting foreground cleaned maps are available 



from the NASA LAM BDA archiveR . Thes e maps are masked 
with the KpO mask jSpergel et al. l2007h that cuts the re- 
gions of brightest galactic emission around the galactic plane 
(~ 20% of the sky), as well as the brightest galactic and 
extragalactic point sources (« 5% of the sky). Zero values 
are assigned to the corresponding pixels. The instrumental 
beam associated with the WMAP radiometers is described 
by an isotropic window function, and the instrumental noise 
is to first order Gaussian, statistically anisotropic, and un- 
correlated. A map with better signal-to-noise ratio can be 
obtained by an optimal combination of the eight foreground 
cleaned and masked maps. At each pixel, this combination 
is obtained by weighting each map by the corresponding in- 
verse noise variance. In order to minimize any error coming 
from the cosmological di pole subtraction, the dipole out- 
side the mask is removed ("Komat su et al.ll2003l '). This over- 
all pre-processing procedure defi nes the so-called t hree-year 
WMAP co-added CMB map pinshaw et al.ll2007h . which is 
used in the subsequent analysis. 

Secondly, ten thousand simulations of the three-year 
WMAP co-added CMB map are considered to compare 
the results of the analysis of the data to what is expected 
from the concordance model. Each simulation is produced 
as follows. Spherical harmonics coefficients of a statistically 
isotropic and Gaussian CMB realization are obtained from 
the angular power spectrum determined by the cosmolog- 
ical parameter s of the three-year WMAP best-fit model 
(jSpergel et al.l l2007') with CAMB (Code for Anisotropics in 
the Microwave Backgrounc^ . The observation at each of the 
eight WMAP radiometers of the Q, V, and W frequencies is 
simulated by convolving that realization in harmonic space 
with the corresponding isotropic window function. Each map 
is then transformed to pixel space at the appropriate resolu- 
tion, and a Gaussian, statistically anisotropic, and uncorre- 
lated noise realization is added with the proper variance per 
pixel. This provides simulations of the CMB, as seen by the 
eight radiometers at the different WMAP frequencies con- 
sidered. The same prescriptions as those described above for 
the data are then applied to produce a three-year WMAP 
co-added CMB map. 

Notice that the WMAP co-added CMB maps for the 
data and s imulations are initi ally produced in HEALPix pix- 
elizatior(f| (|G6rski et al.ll2005h at the resolution N^ide = 512, 
corresponding to maps with more than three million equal- 
area pixels with a spatial resolution of 6.87'. For the sake 
of our analysis, which is applied at 17 scales of the 2GD 
wavelet, corresponding to angular sizes between 27.5' and 
30°, the maps are downgraded to the resolution N^ide ~ 256. 
This provides maps with a bit less than one million equal- 
area pixels with a spatial resolution of 13.7'. Also notice 
that, in pixels close to masked regions, the result of the 
directional correlation of a signal with a steerable wavelet 
is inevitably affected by the zero values of the KpO mask. 
An exclusion mask Ma is therefore deflned at each wavelet 
scale a, identically on the data and simulation s, in order to 
exclude the affected pixels from the analysis (|Vielva et al.l 



^ http://Iambda.gsfc.nasa.gov/ 
^ http://camb.info/ 

http://healpix.jpl.nasa.gov/ 
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Figure 1. Excess kurtosis of the signed-intensity of the three-year 
WMAP co-added CMB data as a function of the 2GD wavelet 
half-width in a range corresponding to angular sizes between 27.5' 
and 30° on the celestial sphere. Data (black rhombi) are com- 
pared with percentiles established from ten thousand statistically 
isotropic and Gaussian simulations produced from the concor- 
dance cosmological model. Significance levels lie roughly below 
1.4% at the four wavelet scales ag, ag, aio, and an, respectively 
corresponding to angular sizes of 8.33°, 10°, 11.7°, and 13.3°. 
The significance level reaches a minimum value of 0.49% at scale 
ag. This identifies a strong detection of non-Gaussianity, in terms 
of an excess of kurtosis in the signed-intensity. 

leaving Na valid pixels from which statistics may be 
estimated. 



3.2 Non-Gaussianity analysis 

The results of the application of the non-Gaussianity analy- 
sis of the local orientation, signed-intensity, and elongation 
to the three-year WMAP co-added CMB map are as follows. 
Let us recall that 17 scales of the 2GD wavelet are probed, 
corresponding to angular sizes on the celestial sphere be- 
tween 27.5' and 30°. The complete list of the wavelet half- 
widths Ohm considered in arcminutes reads: {13.7', 25', 50', 
75', 100', 150', 200', 250', 300', 350', 400', 450', 500', 550', 
600', 750', 900'}. 

Firstly, the significance levels observed for statistics of 
the local orientation and of the local elongation reach mini- 
mum values between 5% and 10%, and only at some isolated 
scales. In other words, these values lie well inside the mid- 
dle region defined above for the distribution of percentiles 
around the median value. We therefore conclude that no de- 
tection is obtained neither for the local orientation nor for 
the local elongation. 

Secondly, the significance levels observed for the excess 
kurtosis of the signed-intensity are roughly below 1.4% at 
the four wavelet scales ag,, ag, aio, and an, respectively cor- 
responding to angular sizes of 8.33°, 10°, 11.7°, and 13.3° 
on the celestial sphere. The excess kurtosis in the signed- 
intensity of the data at each of these wavelet scales is higher 
than the median value defined by the simulations. The sig- 
nificance level reaches a minimum value of 0.49% at scale ag. 
These results identify a strong detection of non-Gaussianity 
in the WMAP co-added CMB data, in terms of an excess of 



Figure 2. Variance of the signed-intensity of the three- year 
WMAP co-added CMB data as a function of the 2GD wavelet 
half-width in a range corresponding to angular sizes between 27.5' 
and 30° on the celestial sphere. Data (black rhombi) are com- 
pared with percentiles established from ten thousand statistically 
isotropic and Gaussian simulations produced from the concor- 
dance cosmological model. Significance levels lie roughly above 
5% between the wavelet scales a2 and 017, corresponding to an- 
gular sizes between 50' and 30° . At the wavelet scale a\ though, 
corresponding to an angular size of 27.5', the significance level 
reaches 0%. This identifies a strong, but isolated detection of 
non-Gaussianity, in terms of a too high variance of the signed- 
intensity. 



kurtosis in the signed-intensity (see Figure [TJ. Also notice 
significance levels of roughly 1.6% and 1.2%, respectively at 
the wavelets scales 04 and 05, corresponding to angular sizes 
of 2.5° and 3.33°. The excess kurtosis in the signed-intensity 
of the data at each of these wavelet scales is well lower than 
the median value defined by the simulations, even though it 
is not considered as anomalous. 

Thirdly, the significance levels observed for the vari- 
ance of the signed-intensity are roughly above 5% between 
the wavelet scales 02 and an, corresponding to angular sizes 
between 50' and 30°, the minimum value being reached at 
scale as. Again, these values lie well inside the middle re- 
gion defined above for the distribution of percentiles around 
the median value. We therefore conclude that no detection 
is obtained at those scales. At the wavelet scale ai though, 
corresponding to an angular size of 27.5', the value of the 
variance is well higher than the median value defined by 
the simulations, and the significance level actually reaches 
0%. Formally, this represents a strong, but isolated detec- 
tion of non-Gaussianity in the WMAP co-added CMB data, 
in terms of a too high variance of the signed-intensity (see 
Figure [2]) . In Section 14] we suggest that it might originate 
in the presence of residual point sources in the data, and 
that it should therefore be discarded. Notice that no detec- 
tion appears neither in the mean nor in the skewness of the 
signed-intensity. 

In summary, the 2GD wavelet gives access to the mea- 
sures of orientation, signed-intensity, and elongation of local 
features of the WMAP temperature data. But a strong de- 
tection of non-Gaussianity is only observed in the excess 
kurtosis of the signed-intensity. This result actually sup- 
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ports the previous detection, with the axisymmetric Mex- 
ican hat wavelet, of an excess of kurtosis in the wavelet co- 
efficient of the WM AP tem perature data (fvielva et al.ll20oi : 
iMukheriee fc Wana.2004 : .Cruz et al-lbOOsf T 



4 SYSTEMATIC EFFECTS 

In this section, we firstly suggest that the high variance of 
the signed-intensity of the WMAP temperature data at the 
smallest wavelet scales, and the corresponding isolated de- 
tection, might originate in the presence of residual point 
sources in the data. Secondly, we discard instrumental noise, 
residual foreground emissions, and large-scale modulations 
of the CMB temperature field related to some unknown sys- 
tematics, as possible origins of the excess of kurtosis in the 
signed-intensity. 



4.1 Residual point sources and lower resolution 

The variance of the signed-intensity of the WMAP co-added 
CMB data is lower than the median value defined by the 
simulations between the wavelet scales og and ai7, corre- 
sponding to angular sizes between 5° and 30° . This variance 
is not considered as anomalous though, as significance levels 
are always roughly above 5%. On the contrary, the variance 
observed between the wavelet scales a\ and as, correspond- 
ing to angular sizes between 27.5' and 3.33°, is higher than 
the median value. Again, the values of the variance between 
the wavelet scales a2 and as are not considered as anoma- 
lous, as significance levels are always roughly above 6%. As 
already emphasized, at the wavelet scale ai , the signifi- 
cance level reaches 0%. The high variance at the smallest 
wavelet scales, in opposition with the behaviour at larger 
wavelet scales (see Figure [5)) , might originate in systematic 
effects such as the presence of residual point sources. Let 
us recall that the three-year WMAP best-fit angular power 
spectrum is obtained after correction for a non-zero best- 
fit amplitude of a re sidual point sources power spectrum 
l|Hinshaw et al.1 120071 '). Consistently, residual point sources 
were recently identified in t he WMAP co-added CMB data 
l|L6pez-Caniego et al.l |2007| ) . These residual point sources 
are however not accounted for in the WMAP co-added CMB 
data analyzed here, while they are accounted for in the sim- 
ulations based on the three-year WMAP best-fit angular 
power spectrum. This could indeed explain the high variance 
observed in the signed-intensity of the data at the smallest 
wavelet scales. Obviously, the contribution of these residual 
point sources is negligible at larger wavelet scales. 

In the absence of a detailed model of the contribution 
of the residual point sources to the variance of the signed- 
intensity, the corresponding isolated detection at the wavelet 
scale ai is discarded, and the analysis is restricted to the 
wavelet scales between ae and an, corresponding to angular 
sizes between 5° and 30° . This does not affect the interpre- 
tation of the detection relative to the excess kurtosis in the 
signed-intensity, which appears between the wavelet scales 
as and an, corresponding to angular sizes between 8.33° 
and 13.3°. For the subsequent analyses dedicated to search 
for the origin of the detection in the excess kurtosis of the 
signed-intensity, the three-year WMAP co-added CMB map 



and the corresponding ten thousand simulations are down- 
graded to the resolution Nside = 32. This provides maps 
with a bit more than twelve thousand equal-area pixels with 
a spatial resolution of 1.83°. For coherence, the initial non- 
Gaussianity analysis is reproduced for the signed-intensity, 
from the data and simulations preliminary downgraded to 
the resolution A^side = 32. All conclusions relative to the first 
four statistical moments remain obviously unchanged. The 
strong detection in the excess kurtosis is slightly enhanced. 
The corresponding significance levels are roughly below 1% 
at the four wavelet scales ag, ag, aio, and an, respectively 
corresponding to angular sizes of 8.33°, 10°, 11.7°, and 13.3° 
on the celestial sphere. The significance level reaches a min- 
imum value of 0.26% at scale ag. 

4.2 Other systematics 

Let us recall the previous detection of an excess of kurtosis in 
the wavelet coefficient of the WMAP tem perature data with 
the axisymmetric Mexican hat wavelet ([Vielva et al.l |2004| : 
[Mukheriec & WanE"2004|; ICruz et aLlbOOSl ). as well as previ- 
ous detections of anomalies obtained in the signed-intensity 
of the WMAP temperature data wit h the 2GD wavelet 
l|McEwen et al.ll2007l : rVielva et al]|2007^ . The corresponding 
analyses concluded that neither instrumental noise nor resid- 
ual foreground emissions are at the origin of the deviations 
observed. These conclusions were drawn from independent 
analyses of the data produced by the eight WMAP radiome- 
ters at the Q, V, and W frequencies. These results suggest 
that the wavelets and statistics used are rather insensitive 
to instrumental noise and residual foreground emissions in 
the WMAP temperature data. In this context, even though 
no similar analysis is performed here, we conclude that nei- 
ther instrumental noise nor residual foreground emissions 
are likely to be at the origin of the excess of kurtosis ob- 
served in the signed-intensity of the WMAP temperature 
data. 

Further investigations are proposed here, considering 
some possible form of unknown systematics. It was recently 
proposed that the WMAP data are possibly affected by a 
large-scale modulation. This modulation was primarily put 
forward as a possible explanation of the North-South asym- 
metry and low mu l tipoles alignm ent of the WMAP data 
('Hel hng et aLll2006l : ICordonllioOTl ). In that framework, the 
WMAP data are of the form T{uj) x [1 -I- /(tj)], where T{uj) 
stands for the CMB temperature on the sky, and ,f{u)) is a 
modulation function containing only low multipoles I. Dipo- 
lar (Z = 1) and dipolar-quadrupolar (Z = {1, 2}) modulation 
functions providing the best-fit cosmological models to the 
three-year WMAP data were propose d. Let us remark tha t 
the best-fit dipolar modulation used (|Eriksen et al.l |2007| ) . 
as well as the best-fit dipolar-quadrupolar modulation used 
ISpergel et al] (2007, see arXiv:astro-ph/0603449vl), were 
not primarily computed for the three-year WMAP co-added 
CMB map itself, but they were shown not to be sensitive to 
the three- year WMAP data set and sky cut. We therefore 
also considered them to be adequate for correction of the 
three-year WMAP co-added CMB map. 

We have checked the stability of the excess of kurtosis 
in the signed-intensity of the three-year WMAP co-added 
CMB data relative to these modulations. Firstly, consider- 
ing the best-fit dipolar-quadrupolar modulation, the strong 
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detection in the excess kurtosis remains unchanged, if it is 
not slightly increased. Just as for the analysis of the non- 
corrected three-year WMAP co-added CMB data, the sig- 
nificance levels are roughly below 1% at the four wavelet 
scales as, ag, aio, and an, respectively corresponding to an- 
gular sizes of 8.33°, 10°, 11.7°, and 13.3° on the celestial 
sphere. The significance levels reach a minimum value of 
0.22% at scales as and ag. Secondly, considering the best-fit 
dipolar modulation, the detection in the excess kurtosis is 
slightly decreased. The significance levels are only roughly 
below 3% at the four wavelet scales ag, ag, aio, and an, with 
a mini mum value of 0.60 % at scale ag. As already empha- 
sized (jVielva et al.ll2007l ). a more precise definition of the 
modulation in terms of specific systematic effects would be 
required before strong conclusions can be drawn from the ap- 
plication of the corresponding corrections. But, even taking 
the results at face value, the proposed dipolar and dipolar- 
quadrupolar corrections are to be rejected as possible origins 
of the observed excess of kurtosis in the signed-intensity of 
the WMAP temperature data. 

In summary, instrumental noise and residual foreground 
emissions are not likely to be at the origin of the excess 
of kurtosis. Large-scale modulations of the CMB related to 
some unknown systematics are explicitly rejected as possible 
origins of the detection. The non-Gaussianity detected in the 
excess kurtosis of the signed-intensity of the WMAP data is 
therefore probably related to the CMB temperature field 
itself. 



5 CONFINEMENT AND DISCUSSION 

We here firstly recall the recent detection of an anomalous 
distribution on the sky of anomalous signed-intensities in 
the three-year WMAP co-added CMB data. We secondly 
test, and tend to reject, the possible confinement of the ob- 
served excess of kurtosis to the directions with an anomalous 
signed-intensity. We finally discuss the detailed interpreta- 
tion of our detections. 

5.1 Local signed-intensity anomalies 

In a very recent analysis of the three-year WMAP co-added 
CMB data with the 2GD wavelet, the distribution on the 
celestial sphere of directions with a signed-intensity anoma- 
lous at 99.865% (formally corresponding to the percentiles at 
three standard deviations {3a) from the mean i n a Gaussian 
distribution) was observed to be anomalous (IVielva et al.| 
I2OO7I ). At the wavelet scale ag, corresponding to an angu- 
lar size of 8.33°, the global significance level of that detec- 
tion, defined as the fraction of the ten thousand simula- 
tions with a number of anomalous directions higher than 
in the data, is 1.39%. The anomalous directions are essen- 
tially distributed in three clusters in the southern galactic 
hemisp here, identify i ng th ree mean preferred directions in 
the skv lVielva et all (|2007l . Figure 5). A first cold spot {i.e. 
with negative signed-intensities) identifies with the anoma- 
lous cold spot originally detected at {9,ip) = (147°, 209°) in 
galactic spheric al coordinates with the axisymmetri c Mexi- 
can hat wavelet (jVielva et al ] |2004l : ICruz et al.|[2005l ) . A sec- 
ond cold spot lies very close to the southern end of the CMB 
dipole axis. The third spot is a hot spot {i.e. with positive 
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Figure 3. Excess kurtosis of the signed-intensity of the three- 
year WMAP co-added CMB data as a function of the 2GD 
wavelet half-width in a range corresponding to angular sizes be- 
tween 5° and 30° on the celestial sphere. Data (red squares) are 
compared with percentiles established from ten thousand statisti- 
cally isotropic and Gaussian simulations produced from the con- 
cordance cosmological model. The directions with an anomalous 
signed-intensity are excluded from the statistical analysis, in the 
data as well as in each of the ten thousand simulations inde- 
pendently. Significance levels are roughly below 1% at the four 
wavelet scales ay, ag, ag, and aio, respectively corresponding to 
angular sizes of 6.67°, 8.33°, 10°, and 11.7°. The significance level 
reaches a minimum value of 0.28% at scale ag. This still identifies 
a strong detection of non-Gaussianity, in terms of an excess of 
kurtosis in the signed-intensity. 

signed-intensities) close to the southern end of the eclip- 
tic poles axis. The detection is confirmed at the neighbour 
wavelet scales ag, aio, and an, respectively corresponding 
to angular sizes of 10°, 11.7°, and 13.3°. Instrumental noise, 
residual foreground emissions, as well as large-scale modu- 
lations of the CMB related to some unknown systematics, 
are rejected as possible origins of the detection. The local- 
ized anomalous distribution of anomalous signed-intensities 
identified may therefore probably be imputed to the CMB 
temperature field itself. 

5.2 Confinement analysis 

Postulating the Gaussianity of the CMB, one may interpret 
the anomalous distribution of directions with an anomalous 
signed-intensity as a clear departure from statistical isotropy 
(|Vielvaet al.ll2007l ). But more generally, the anomaly ob- 
served highlights a deviation of the CMB temperature field 
from the whole assumption of statistical isotropy and Gaus- 
sianity. In the context of the present non-Gaussianity de- 
tection, the anomalous distribution of anomalous signed- 
intensities previously identified represents a serious candi- 
date to explain the excess of kurtosis observed in the signed- 
intensity of the three- year WMAP co-added CMB data. This 
idea is supported by the fact that the two detections are ob- 
served at the same wavelet scales. The hypothesis to check 
is to know if the whole non-Gaussianity observed is confined 
to the directions with an anomalous signed-intensity, in the 
idea that the excess in the number of anomalous directions 
would bear the departure from Gaussianity. 
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The confinement analysis simply consists in reproduc- 
ing the previous analysis on the first four statistical moments 
of the signed-intensity of the three-year WMAP co-added 
CMB data, from which the directions with an anomalous 
signed-intensity are excluded. The only way to account for 
a possible bias introduced by the exclusion of the extremal 
values of the data above a given threshold is to apply the 
same exclusion process to each simulation independently. 
In particular, the coherence of the procedure can only be 
achieved if all values above the threshold, and not only part 
of it, are identified and excluded from the statistical analysis, 
in the data and in the simulation^. Notice that, while the 
anomalous signed-intensities were originally identified at a 
threshold of 99.865%, we consider here a threshold at 99.5%. 
This lowering is performed in order to avoid a possible neg- 
ative conclusion, relative to the confinement, simply due to 
the fact that directions in the immediate vicinity of the di- 
rections thresholded are not excluded from the analysis. For 
completeness, even though the detections are only observed 
at the four wavelet scales between ag and an, corresponding 
to angular sizes between 8.33° and 13.3°, the confinement 
analysis is performed at all wavelet scales between ae and 
ai7, corresponding to angular sizes between 5° and 30°. 

The results of the analysis performed are as follows. 
The strong detection of an excess of kurtosis in the signed- 
intensity of the WMAP temperature data is preserved. The 
significance levels are roughly below 1% at the four wavelet 
scales 07, as, ag, and aio, respectively corresponding to an- 
gular sizes of 6.67°, 8.33°, 10°, and 11.7° on the celestial 
sphere. The significance level reaches a minimum value of 
0.28% at scale as (see Figure |3]). The variance of the signed- 
intensity at each wavelet scale was well lower than the me- 
dian value defined by the simulations before the threshold- 
ing. It was not considered as anomalous as the significance 
level reached a minimum value of roughly 5% at scale ag. 
The variance of the signed-intensity at each wavelet scale is 
still well lower than the median value defined by the sim- 
ulations after the thresholding. It is still not considered as 
anomalous, even though the significance level reaches a min- 
imum value of roughly 1% at scale ag. No detection ap- 
pears neither in the mean nor in the skewness of the signed- 
intensity. 

In summary, removing the directions with an anomalous 
signed-intensity does not solve the observed discrepancy be- 
tween the WMAP temperature data and simulations. Con- 
sequently, taking the CMB temperature angular power spec- 
trum of the concordance cosmological model at face value, 
we can conclude that the strong detection in the excess kur- 
tosis of the signed-intensity of the WMAP temperature data 
is not confined to the directions with an anomalous signed- 
intensity. 



^ We therefore have to assume that, in the data, all extremal 
values of the background Gaussian distribution above the thresh- 
old are excluded. The only formal reason for which extremal 
values outside the effectively observed directions with a signed- 
intensity above the threshold might have been missed would be 
that some non-Gaussianity coincidentally compensates for the ex- 
tremal value in the corresponding direction on the celestial sphere. 
In that trivial case, one readily knows that the non-Gaussianity 
observed in the WMAP data is not confined to the localized dis- 
tribution of anomalous signed-intensities. 



5.3 Discussion 

Firstly, if the excess of kurtosis observed in the signed- 
intensity of the WMAP temperature data is related to 
the CMB temperature field itself (see Section |4} , the fact 
that this non-Gaussianity is not confined to the direc- 
tions with an anomalous signed-intensity seems natural. In- 
deed, independently of the modification of the basic infla- 
tionary scenario that might exp lain the non-Gauss ianity 
of the CMB temperature field (jBartolo et all |2004 ), the 
cosmological principle still implies its statistical isotropy, 
at least as a first approximation. Non-Gaussian pertur- 
bations are therefore more naturally widely spread over 
the whole sky. We also notice that the wavelet scales at 
which the non-Gaussianity is observed are compatible with 
the size of CMB anisot ropics due to topolog ical defects 
such as cosmic textures (|Turok fc Spergellll990l ). or due to 
the Integrated Sachs-Wolfe effect, which is associated with 
the time evoluti on of the gravitational potential of large 
scale structures JSachs fc Wolfelll967l : iRees fc Sciamalll9"68l : 
iMartmez-Gonzalez fc Sand Il990l 'l. Texture models suggest 
the presence of a number of textures with angular sizes above 
1°, which can induce hot s pots or cold sp ots of corresponding 
angular size in the CMB llTurok fc Spergel 19901. A recent 
Bayesian analysis (ICruz et al.l 2007 ) showed that the cold 
spot originally detected at {9,ip) = (147°, 209°) in galactic 
spherical coordinates with the axisymmetric Mexican hat 
wavelet, is satisfactorily described by a texture with an an- 
gular size on the celestial sphere around 10° . Other analyses 
also showed that the time evolution of the gravitational po- 
tential of large scale structures such as voids might induce 
cold spots in the CMB with angular sizes of several degree s 
on the celestial sphere (|Martmez-Gonzalez fc Sand Il990l ) . 
The cold spot identified at {6,ip) = (147°, 209°) in galactic 
spherical coordinates with an angular size around 10° could 
actually be explained in terms of a void at a redshift 2 ~ 1 
and with a diameter arouri d 300/i"^ Mpc (|Inoue fc Silk.20061 , 
120071 : iRudnick et al.|[2007h . 

Secondly, as already emphasized, an excess of kurtosis 
in the wavelet coefficient of the WMAP temperature data 
was previously detected with the axisymmetric Mexican hat 
wavelet. The non- Gaussian deviation observed with the ax- 
isymmetric Mexican hat wavelet is undoubtedly related to 
the present detection in the excess kurtosis of the signed- 
intensity with the 2GD wavelet, notably because it is ob- 
served with a similar statistics at the same angular sizes 
on the celestial sphere. From this point of view, both de- 
tections support one another. But the axisymmetric Mexi- 
can hat wavelet also allowed the detection of the cold spot 
at {9,ip) = (147°, 209°) in galactic spherical coordinates, 
which was interpreted t o be the exclusiv e origin of the ex- 
cess of kurtosis detected l|Cruz et al.ll2005l ') . On the contrary, 
we have concluded that the detection in the excess kurto- 
sis of the signed-intensity of the WMAP temperature data 
with the 2GD wavelet is not confined to the previously iden- 
tified directions with an anomalous signed-intensity. Con- 
sequently, even though the two detections are similar, they 
probably simply do not identify the same non-Gaussian con- 
tent in the WMAP temperature data. 

Finally, let us also underline that the values of the cos- 
mological parameters are affected by uncertainties associ- 
ated with the limited precision of measurement of the CMB 
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temperature angular power spectrum. These uncertainties 
are associated with the cosmic variance, but also with sys- 
tematic effects such as instrumental noise and residual fore- 
ground emissions. The simulations produced for our anal- 
ysis are obtained from the angular power spectrum deter- 
mined by the cosmological parameters of the concordance 
model (i.e. the three- year WMAP best-fit model). They do 
not account for the quoted uncertainties. Consequently, be- 
fore giving full credit to our conclusions, a deep analysis 
should be performed to check the stability of the various 
detections considered when the WMAP temperature data 
are compared with simulations produced from any possi- 
ble angular power spectrum inside the experimental error 
bars. Formally, any of our conclusions might be affected by 
this further analysis, from the fact that the non-Gaussianity 
observed in the WMAP temperature data is not confined 
to the directions with an anomalous signed-intensity, up to 
the mere detection of an excess of kurtosis in the signed- 
intensity. However, such an analysis would be very involved 
and is not produced here. On the one hand, the stability of 
the detection of an excess of kurtosis in the wavelet coeffi- 
cient of the WMAP temperature data with the axisymmet - 
ric Mexican hat wavelet was sugg ested (|Vielva et al.ll2004 ). 
The same conclusion probably holds for the present detec- 
tion of an excess of kurtosis in the signed-intensity of the 
WMAP temperature data with 2GD wavelet. On the other 
hand, the confinement analysis itself is based on the pre- 
vious detection of the dist ribution of direct ions with an 
anomalous signed-intensity l|Vielva et al.ll2007h . At present, 
no analysis confirmed the stability of this distribution rel- 
ative to the uncertainties on the cosmological parameters. 
A possible excess of power in the concordance model rela- 



tive to the WMAP t emperature data ( Spergel et al.l l2003l : 
iHinshaw et all l2007l : iMonteserin et al.1 120071 ') might imply 



that a part of the distribution on the celestial sphere of 
directions with an anomalous signed-intensity was actually 
not detected. This would probably not question the fact that 
this distribution is anomalous at the wavelet scales between 
as and an, corresponding to angular sizes between 8.33° 
and 13.3°, but would simply suggest that the global signif- 
icance level for the detection was underestimated. Ifowever 
in such a case, the confinement analysis itself, which explic- 
itly requires the exclusion of all the extremal values above 
a given threshold, both in the data and in the simulations, 
might not be performed anymore. No conclusion relative to 
the possible confinement of the non-Gaussianity observed 
to the directions with an anomalous signed-intensity could 
therefore be reached. 



by comparison of the data with ten thousand statistically 
isotropic and Gaussian simulations produced from the con- 
cordance cosmological model. No detection is made neither 
in the orientation analysis nor in the elongation analysis. A 
strong detection is made in the excess kurtosis of the signed- 
intensity of the WMAP data, with a significance level below 
0.5% at a wavelet scale corresponding to an angular size 
around 10°, and confirmed at neighbour scales. This sup- 
ports a previous detection of an excess of kurtosis in the 
wavelet coefficient of the WMAP data with the axisymmet- 
ric Mexican hat wavelet. An isolated detection is also made 
in the variance of the signed-intensity at the smallest wavelet 
scale. Systematic effects such as residual point sources in the 
WMAP co-added CMB data are suggested to originate this 
anomaly, which is consequently simply discarded. 

Instrumental noise and residual foreground emissions 
are not likely to be at the origin of the detection in the 
excess kurtosis of the signed-intensity. Large-scale modu- 
lations of the CMB related to some unknown systematics 
are explicitly rejected as possible origins of the detection. 
The observed non-Gaussianity is therefore probably to be 
imputed to the CMB temperature field itself, thereby ques- 
tioning the basic infiationary scenario upon which the con- 
cordance cosmological model relies. In this context, taking 
the CMB temperature angular power spectrum of the con- 
cordance cosmological model at face value, further analy- 
sis also naturally suggests that this non-Gaussianity of the 
WMAP temperature data is not confined to the localized 
distribution of anomalous signed-intensities. However, this 
last result, in particular, might be sensitive to uncertain- 
ties affecting the cosmological parameters. Further analyses 
should be performed before giving it full credit. 
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6 CONCLUSION 

The decomposition of a signal on the sphere with the steer- 
able wavelet constructed from the second Gaussian deriva- 
tive gives access to morphological measures such as the ori- 
entation, signed-intensity, and elongation of the signal's lo- 
cal features. In this work, the three-year WMAP co-added 
data of the CMB temperature field are analyzed through 
the first four statistical moments of the random fields asso- 
ciated with these local morphological measures, at wavelet 
scales corresponding to angular sizes between 27.5' and 30° 
on the celestial sphere. The statistical analysis is performed 
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